Active Width at a Slanted Active Boundary in Directed Percolation 
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The width W of the active region around an active moving wall in a directed percolation process 
diverges at the percolation threshold p c as W ~ At~"W ln(^), with e = p c — p, eo a constant, and 
vn = 1.734 the critical exponent of the characteristic time needed to reach the stationary state 
£ii ~ e -I/ H . The logarithmic factor arises from screening of statistically independent needle shaped 
sub clusters in the active region. Numerical data confirm this scaling behaviour. 



I. INTRODUCTION 

Directed percolation (DP) has emerged as one of the 
generic absorbing state type dynamic processes. It de- 
scribes epidemic processes, e.g., forest fires and various 
types of surface catalysis processes (l]-|] . Such processes 
include a so-called absorbing state, typically the vacuum, 
from which it can not escape. The relevant tunable pa- 
rameter is the propagation probability p. The system un- 
dergoes a phase transition from the absorbing phase at 
small p, where the stationary state is the absorbing state, 
into an active stationary phase at large p, where the sys- 
tem refuses to die. The scaling properties at DP dynamic 
phase transitions are known for almost two decades, and 
it's now realized that DP critical behaviour is the generic 
universality class for dynamic absorbing state type pro- 
cesses 0. 

At DP type critical points the equilibration time £|| 
diverges. It scales as £|| ~ £J compared to the spatial 
correlation length with dynamic exponent z = 1.581 
||. For example, starting from a single seed, the survival 
probability obeys the scaling form 



P s (e,t) = b-*°P s (b 1 ^e,b- z t) 



(1) 



with e = p c — p the distance from the critical point. This 
leads to 



-t . 



P s ~e"exp(— ), 



(2) 



with exponent /3 = x s vj_. The exponential factor reflects 
that deep inside the absorbing phase P s decays exponen- 
tially in time. The equilibration time diverges at the DP 
critical point as £|| ~ e~ v w with z = vn/v±. At p c the 
survival probability decays as a powerlaw, P s (t) ~ t~ s 
with S — x s /z = /3/vu. 

A recent direction of research in this topic concerns the 
scaling properties near boundaries - 10 1 . Those studies 
address absorbing and reflective walls. The scaling prop- 
erties are modified by surface type critical exponents. In 
particular, the survival probability for a seed near the 
boundary obeys the same scaling form as above, but with 



a new interface critical exponent x, and therefore a mod- 
ified value for (3. 

In this study we discuss the scaling properties near 
active boundaries. Consider a stationary active vertical 
wall in the system. All sites in the wall are alive. The 
critical exponent j3 is not an issue, because the system 
remains active near the wall for all p. However, in the 
absorbing phase the cloud of active sites near the wall 
has a specific stationary state width, which is expected 
to diverge as W ~ £j_ ~ e~ v± . Widths like this diverge 
with bulk exponents. 

Assume that this wall is slanted, with an arbitrary an- 
gle 6 ^ 90° with respect to the horizontal direction (see 
Fig.|l|). In the space-time interpretation of the configu- 
rations, the wall moves with a constant velocity. It acts 
as a slanted active curtain rod. A curtain of active sites 
hangs down from it as illustrated in Fig.|l|. For p < p c the 
curtain has a finite width l± and length {m — l± tan(#). 

In this study we address how the stationary state width 
of this slanted curtain scales near the DP critical point. 
Naively this seems a simple question. One would expect 
that the curtain width diverges with the same exponent 
as the equilibration time scale, W ~ e _I/ ii , i.e., with the 
same exponent as the length of a curtain hanging down 
from an horizontal curtain rod (6 = 0). The latter is 
equivalent to asking for the survival probability in the 
set-up without any walls where all sites are active in the 
initial state. 

This expectation is based on the anisotropic scaling 
properties. Consider a system with a rod at angle 9 =/= 0. 
The horizontal and vertical bulk lengths diverge with 
different exponents, as £|| ~ Therefore, a system 

at p c — p = e and wall angle 9 is equivalent by renor- 
malization to a system with a smaller wall angle 6' at 
e' = fc-V^e with tan(0') ~ b^ 1 tan(0). The scaling 
properties of W should not depend on the angle 9, since 
the rod renormalizes towards the horizontal position. We 
should expect the same scaling behaviour as at 9 = 0. 
However, a recent numerical study |i"lf| seems to contra- 
dict this. 

Kwon et. al. [ pd| studied a model with two absorbing 
states. It undergoes a dynamic phase transition which 
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belongs to Directed Ising (DI) universality class when 
the two absorbing states are symmetric, and belongs to 
the Directed Percolation (DP) universality class when 
a symmetry breaking field is introduced. They stud- 
ied the interface dynamics of the active domain between 
two asymmetric absorbing states. As one absorbing state 
dominates over the other, the interface is driven into the 
unpreferred absorbing region with a constant velocity. 
Therefore they expected the width of the active domain 
to scale like the horizontal width of the active curtain 
in the above setup for ordinary DP models. A simple 
power-law fit of their data suggests that the active do- 
main width scales as W ~ e~ x with x ~ 2.00(5), which 
does not agree with the DP exponent v\\ ~ 1.734. 

In this paper we address the same issue more directly. 
We insert a slanted active wall into the most basic model 
for DP, the one studied orig inally by Kinzel (HH, see 
section 2. We find a similar anomalous value for the 
width exponent. W ~ e~ x scales as x ~ 1.95(5). In 
section 3 we develop a qualitative scaling theory. It pre- 
dicts that the curtain width scales with the conventional 
exponent v\\ but with an additional logarithmic factor as 
W ~ At~ v w ln(^). In section 4 we show that the numer- 
ical Monte Carlo data fits this form well. In section 5 
we illustrates how DP type processes with slanted walls 
can be studied in the Master equation formalism. Our 
finite size scaling (FSS) results, using exact numerical 
enumeration of the eigenvalue spectrum, show that at p c 
the width of the slanted curtain diverges as W ~ L z with 
system size. This confirms the absence of a new indepen- 
dent exponent. The logarithmic factor arises only in the 
e dependence. 

II. NUMERICAL RESULTS FOR THE CURTAIN 
WIDTH. 

Consider the square space-time lattice shown in Fig- 
ure |. All bonds run under 45 degrees. The black (open) 
circles represent the active (inactive) sites. Time evolves 
from top to bottom in half units t — > t + \ . Bonds be- 
tween nearest neighbor sites at t and t — » t + 4 are being 
created with probability p but only if the upper site is 
active. Each bond activates the lower site. Kinzel stud- 
ied this model in detail with Master equation type FSS 
in the early eighties [jl2|. The critical exponents and the 
location of the DP transition are known quite accurately. 
For example, the latest series expansion results put the 
DP phase transition at p c w 0.6447 

We modify the boundary conditions in this model to 
accommodate an active wall. The lattice is semi-infinite, 
bound to the left by the wall, which runs away under 
6 = 45° as shown in Fig. ^. 45° is its natural angle 
for the curtain rod for this specific lattice. We can re- 
strict ourselves to this angle because the scaling proper- 
ties of the curtain width should not depend on the angle 



according to the anisotropic scaling argument outlined 
above. Moreover, the angle is a continuous parameter in 
the model by Kwon et. al. |ll| and their results show no 
angle dependence. 

We perform Monte Carlo simulations with as initial 
configuration an active wall in an inactive bulk. The hor- 
izontal curtain width is defined as the distance of the last 
active site from the rod in each time slice. For p < p c , the 
width grows initially approximately linear in time, until 
it saturates at the stationary state value which varies 
with e = p c — p. Figure |3| shows the active width versus 
e on a logarithmic scale. The line is quite linear over the 
two decades shown. The slope is clearly distinct from the 
expected value vn sa 1.734 and close to the value found 
by Kwon et. al. jo). In Fig.^ we perform a more care- 
ful FSS analysis to the same data. We fit the numerical 
data from two nearby points, ti — y/2ei, to the form 
W ~ ae~ x and plot function of e, the exponent 

x appears to be around 1.95. This fit is remarkably sta- 
ble, and shows virtually no power law type corrections to 
scaling. Taken out of context it is strongly suggestive of 
a new independent critical exponent. The other curves in 
Fig.|| relate to the FSS analysis assuming an additional 
logarithmic factor as discussed in the next two sections. 



III. INDEPENDENT CLUSTER 
APPROXIMATION 

Figure [l] shows a typical curtain configuration in a 
Monte Carlo simulation at a p just below the percola- 
tion threshold p c . The most striking features are the 
needles in the curtain. Isolated clusters are expected to 
be needle like. The correlation length in the time di- 
rection diverges faster than in the spatial direction, as 
£y ~ Therefore, active clusters (when grown from a 
single seed) become needle shaped near the percolation 
threshold. Figure gives the impression that close to p c , 
the curtain consists of a set of weakly interacting needle 
shaped clusters when viewed from length scales larger 
than £x- 

In this section we pursue the implications of the as- 
sumption that such needles are completely uncorrelated. 
In that approximation the probability that the curtain 
extends over a horizontal distance I is given by the prob- 
ability that a needle longer than r = I tan(0) hangs down 
from the curtain rod vertically above that site. Let P be 
that probability. It must have the same form as the sur- 
vival probability from a single seed, Eq. (gj) , The actual 
value of the exponent (3 turns out to be irrelevant in this 
section, but it must be identical to the single seed value, 
according to a time reversal symmetry argument [ [l4| . 

The spatial coordinate needs to be coarse grained, 
because the needles can only be uncorrelated beyond 
the horizontal correlation length £j_ ~ e~ v± . Define 
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n = x/£± as the coarse grained discrete spatial coor- 
dinate and recall that t = xtan(#) is the corresponding 
vertical distance from the curtain rod to the same point. 
The probability for the curtain to have width n factorizes 
in the independent needle approximation as 



P w (n) = P(n) ] ; 



1 - P(n') 



(3) 



This equation can be rewritten into a derivative form 
P w (n + 1) - P w (n) P(n + 1) - P(n)[l - P[n + 1)] 



P w {n + 1) 



P(n + 1) 



(4) 



The maximum of the distribution obeys the relation 

P w {n-l) = P w {h) (5) 

and can be written as 

1 1 



P(n) P(h-1) 



1. 



(6) 



Assume that P has the same asymptotic form as the sin- 
gle seed survival probability, in Eq. (J2J) , and that the max- 
imum of the distribution occurs in this range of n. The 
transformation to the coarse-grained n = x/£± ~ xe u± 
variable changes the critical exponent inside the expo- 
nential factor 



P ~ Be (i e- hnt 



(7) 



with A = V\\ — v±, and b ~ tan(0). Inserting this form 
into Eq.(0) leads to 



1 



(8) 



and, after expanding the exponential on the left hand 
side, to 



6n|e| A ~ln(^) + (/?-A)ln(e) 
o 



In original units this reads 



(9) 



(10) 



The characteristic probability depends on the wall an- 
gle as eq ~ 1/ tan(0). The most probable width W scales 
with the expected exponent i^y but contains an additional 
logarithmic factor. 

Asymptotically the most probable and the average 
widths coincide. Eq.(|]) can be approximated in the con- 
tinuum limit as 



1 dP w 

P,„ dn 



= 1 - 



P(n) 

P(n + 1; 



+ P{n) 



(11) 



Close to p c and for large n, where P obeys Eq.(^), we 
can integrate this 



P„(n) 



exp[(l 



e»^)n 



B 



B 



e 0-A e -bne^ 



(12) 



expL--e^- A e- 6n ^ 



This distribution decays exponentially on both sides of 
the most probable value and becomes sharp at the crit- 
ical point, e — > 0. We checked explicitly that the most 
probable and average coincide in this limit, and scale 
asymptotically with the same logarithmic factor, as in 
Eq.©. 



IV. LOGARITHMIC CORRECTIONS TO 
SCALING ANALYSIS. 

The logarithmic factor in the independent needle ap- 
proximation formula for the curtain width 



W(e) ~ Ac'"" ln(-) 
e 



(13) 



does not change the asymptotic exponent. It is still equal 
to mi . However the finite size scaling (FSS) approach to 
this value is very singular. A conventional FSS analysis 
involves the construction of approximants for the critical 
exponent x by fitting the values of W at to nearby e to 
a pure power law form, W ~ e _x . This is equivalent 
to defining x(e) as a derivative and yields for the above 
logarithmic form 



e dW _ 1 
X ~ ~W^k ~ U]l + mpO 



(14) 



This function approaches f|| in a singular manner. In the 
interval 0.01 < e/e < 0.3, x seems to converge convinc- 
ingly with a linear correction to scaling term to an effec- 
tive exponent which is about 0.2 too large. One would 
have to go to extremely small e's to see the true conver- 
gence. The power law fit in Figfl shows signs of this. 

The two other curves in FigMshow the FSS estimates 
for the exponent v\\ according to the form Eq.(|l3|) with 
6q = 1 or eq = 0.5. eo is unknown, but likely of or- 
der one. Both curves converge towards the conventional 
value v\\ = 1.734. This is strong evidence for the presence 
of the logarithmic factor. 



V. FINITE SIZE SCALING AT THE 
PERCOLATION THRESHOLD 

The logarithmic factor originates from the screening 
of independent needles. It should not play a role in the 
FSS at the percolation threshold itself, because there £j_ 
diverges, and the independent needle concept becomes 
meaningless. 
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So the curtain width must scale as W ~ L z at p c , if it is 
really true that no independent new exponent is involved. 
To confirm this we present in this section numerical data 
from master equation type finite size scaling using ex- 
act enumeration. We also performed Monte Carlo sim- 
ulations but prefer to present our master equation data 
since this method requires a technical novelty. 

A moving wall is inconvenient in simulations. The lat- 
tice is finite by necessity and the moving wall requires 
a much bigger lattice than the one actually used by the 
process. This is a handicap in particular for master equa- 
tion calculations where one evaluates the rate at which 
the stationary state is being reached by letting time go 
to infinity at each lattice size L. Those systems sizes are 
typically small, L < 20 in our case, because phase space 
scales exponentially with L. Compared to MC simula- 
tions the master equation method trades system size for 
numerical accuracy, and the ability to perform a detailed 
corrections to scaling analysis. The accuracy of the two 
methods is typically comparable, except for specific is- 
sues, like the logarithmic factor in the previous sections, 
which require intrinsic large lattice sizes. 

The solution to the moving wall problem is to distin- 
guish between the time and spatial directions of the dy- 
namic process, e± and en, and the ones used in the mas- 
ter equation. There is no need for them to coincide. We 
choose a set-up where the master equation's time and 
space directions arc redirected in the following manner. 
Lines of constant time are parallel to en — ej_, such that 
the moving wall coincides with the t = line. Lines of 
constant position are parallel to the x-axis, which in the 
dynamic process represented lines of constant time. 

The following skewed dynamic rule implements this 
pace-time rotation. Consider a square space time lattice 
(Fig.|| rotated over 45 degrees). Each site in the master 
equation time slice t is updated sequentially from right 
to left. The probability for site x at time r to be active 
depends on whether site x — 1 was active at the previ- 
ous time t — 1 and/or at this moment in time, t. This 
set-up requires screwed boundary conditions. The forest 
fire runs under an angle. In this new interpretation the 
active wall represents a fully active initial configuration. 

The energy gap in the spectrum of the time evolution 
operator (transfer matrix), is related to the curtain width 
in the following manner. Let \I) be the initial state of 
the master equation, |0) the absorbing state, and T be 
the transfer matrix. The stochastic nature of the trans- 
fer matrix implies that the disordered state \D) is a left 
eigenvector with eigenvalue Ao = 1. Define (Xx the 
projection operator which returns one (zero) when site 
x active (inactive) . The curtain width is associated with 
the probability distribution for site x to be active at time 
t but after that never again. This takes the form 

P(x,t)= lim (D\ [(1 - a x )f a x f *\I) (15) 

tp — >oc 



The operator (1 — a x )T has Ao = 1 as largest eigenvalue 
since 

(l-a x )f]\Q) = \Q) (16) 

and because attaching a projection operator to T can not 
result in an eigenvalue larger than the largest in T. Let 
(L x | be the corresponding left eigenvector (which can be 
evaluated numerically). Inserting this leads to 

P{x,t) = (L x |a x f*|/) (17) 
~ (£ a |a*|Ai)At<Ai|i) 
- exp[-t/f t ] 

with £t — log(Ai) and Ai the next largest eigenvalue of 
f. 

This illustrates that the curtain width scales in the 
same manner as the the characteristic time £t needed to 
reach the stationary state, when the latter is measured in 
this space-time twisted coordinate system. Fig.[s] shows 
the FSS estimates for the dynamic exponent z according 
to £ t ~ L z and W ~ L z . Both converge clearly to the 
DP dynamic exponent z = 1.58. This confirms that no 
new independent curtain width exponent is presents. 

VI. FINAL REMARKS 

The analysis presented in this paper explains the 
anomalous scaling of the width of the slanted curtain 
boundary in DP type processes. The needles screen each 
other, and that leads an extra logarithmic factor accord- 
ing to the independent needle approximation. Our nu- 
merical data confirm the validity of this assumption. 

The same mechanism must apply to other dynamic 
processes, like directed Ising type absorbing state dy- 
namics, and also to other quantities. Consider the fol- 
lowing example. Directed percolation describes epidemic 
growth processes without immunization, where the prob- 
ability to be sick at time t + 1 requires that you yourself 
or at least one of your neighbours is already sick at time 
t. Consider an initial condition that everybody is sick at 
time t = 0. A stationary local observer will conclude that 
below the percolation threshold the life time of the epi- 
demic scales as t ~ e _!/ ll . A moving observer concludes 
it diverges faster, as t ~ e~ v w log(^). 
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FIG. 2. Lattice structure near the active boundary 




FIG. 1. The curtain of active sites at the active slanted 
boundary 
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FIG. 3. Log plot of active width versus p c —p from straight 
Monte Carlo simulations on unlimited system sizes. The solid 
line represents the data. The dashed straight lines of slopes 
— 2 and —1.734 are guides to the eyes. 
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FIG. 4. Estimates for the active width exponent, x. In fit 
(a), W is assumed to scale as W ~ e~ x , in (b), as W ~ e~ x In e 
and in (c), as W ~ e~ x (lne + In 2). 
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FIG. 5. Finite size scaling exponent z for (x) the charac- 
teristic active width, W ~ L z , and for (+) the time to reach 
the stationary state, t ~ L z , at the percolation threshold in 
the transfer matrix set-up. The data virtually coincide 
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